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Abstract 

In view of their physics goals, future neutrino factories from muon decay aim 
at an overall flux precision of 0(1%) or better. We analytically study the QED 
radiative corrections to the neutrino differential distributions from muon decay. 
Kinematic uncertainties due to the divergence of the muon beam are considered 
as well. The resulting corrections to the neutrino flux turn out to be of order 
0(0.1%), safely below the required precision. 
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1 Introduction 



Results on neutrino oscillations from Superkamiokande pQ and SNO |2] provide a com- 
pelling evidence for neutrino masses, constituting the first strong indication of Physics 
beyond the Standard Model. Much is still unknown, though, regarding fundamental 
issues such as the absolute neutrino mass scale, the possible Majorana character of 
neutrino fields, the ordering of their mass eigenstates with respect to charged lepton 
eigenstates, or the possible existence of leptonic CP violation and its tantalizing rela- 
tionship to baryogenesis. In this situation one could argue that the subject of lepton 
flavor physics is at its exciting infancy, and to obtain rough answers to those questions 
could be a sufficient goal at present, postponing any aim at a precise determination of 
the involved parameters. Nevertheless, some of those questions prerequire precision: 
for instance the study of CP violation rests upon a precise knowledge of the angles in 
the neutrino mixing matrix. 

In a more general way and much as for the quark sector, it is necessary to know 
accurately the values of the masses and mixing parameters in the lepton sector, as a 
first step to unravel the flavor puzzle. And what does precision means, quantitatively?. 
For instance, with which precision is it desirable to determine the values of the leptonic 
mixing angles in order to discriminate between models for neutrino masses? Clearly 
no definite answer can be given to such question, but as an indication it has been 
argued jl] that a 10% — 1% precision in the knowledge of, say, sin 2 29 atm would result 
in significant advance 1 . It is not impossible to envisage such a precision. In resume, we 
are simultaneously entering a discovery and a precision era in neutrino physics. With 
the bonus that the extraction of physical conclusions will not be necessarily hindered 
by large theoretical errors, as it happens in the quark sector due to QCD long distance 
contributions. 

A quest for precise physics answers evidently requires an effort in precision on 
the experimental conditions, and on the knowledge of the neutrino flux to start with. 
Several experiments using neutrino beams from particle accelerators such as K2K, 
MINOS and OPERA 3J will take data in the next few years. Their reach will be limited 
by the use of conventional neutrino beams produced from a charged pion source. The 
decay n + — > /U + z/ M (n~ — > n^u^) produces a beam with a 0(1%) component of v e 
from kaon decays. The v e contamination limits the precision of the flux measurements, 
resulting in an error of 7% for K2K, while MINOS reduces it to 2% A further 
step forward could be provided by the so-called superbeams which, although based on 
the same traditional beams, can achieve better precision thanks to the much higher 
statistics. It has been argued, for example, that by working at energies below the 
threshold of kaon production, the v e flavor contamination could be reduced, with the 
overall figure of merit for precision in the flux measurements limited to (9(1%) [H1IZ1- 

1 9atm denotes the mixing angle dominantly responsible for the atmospheric oscillations, denoted 
by #23 in the by now standard parameterization |S] 
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A major advance should come from a neutrino factory from muon decays, aiming 
at both fundamental discoveries and (9(1%) precision measurements. Present projects 
consider the production of very intense muon sources of about 10 20 muons per year jH] . 
Neutrino beams originate from the decay of high-momentum muons along the straight 
sections of a storage ring. The beam produced presents a precisely known neutrino 
content: 50% muon neutrinos and 50% electron antineutrinos if a \i~ beam is used, 
and 50% muon antineutrinos and 50% electron neutrinos if a /i + beam is used. The 
resulting v fluxes are expected to be known with a precision better than 1% Pj. It is 
necessary to ensure that any possible corrections and sources of errors are controlled 
at that level. In this work, we study two effects: the contribution of QED one-loop 
corrections to muon decay and the divergence of the muon beam. For both cases, we 
give novel corrected formulae for neutrino differential distributions. 

Radiative corrections to the electron differential distribution in fj,~ — > e~ + v e + 
were calculated long ago resulting in a correction of 0(1%) jTU], larger than the 
expected effect of O(-) ~ 0(0.1%). Such an effect is at the level of the expected 
precision at a neutrino factory. In this work we study whether QED corrections affect 
neutrino distributions at the same order. 

The correction to the (massive) neutrino spectra from unpolarized muons has been 
first calculated in jTTj. In our work, we give new anlytic formulae with m u = and 
m e = including muon polarization, relevant for neutrino factory measurements. Dif- 
ferent from the electron case, the analysis of the correction to the neutrino differential 
distributions entails non-trivial integrations. By using the correspondence between the 
QED corrections to the /z-decay and those for the charge 2/3 heavy quarks in QCD, 
we make use of the techniques developped in Refs. O HS1 El for the calculation of 
the QCD corrections to the lepton spectrum in the decay t — > b + / + + v\. 

The second subject addressed in this paper is that of the muon beam divergence, 
one of the basic properties that can bias the predicted neutrino spectra . We explore 
the error induced in the neutrino distributions at the far site due to the systematic 
uncertainty on the angular divergence, and compare our results with previous ones in 
which this effect was not included |15j . 

The paper is organized as follows. In section 2 we recall the tree-level angular 
distributions. In section 3 the neutrino one- loop corrected formulae are given, with 3.1 
specializing in the soft photon limit and cancellation of infrared divergences. Section 
4 accounts for the corrections due to the beam divergence. 
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2 General definitions 



In the muon rest-frame, the angular distributions of the neutrinos produced in the 
decay — > e~ + + u e , Fig. [T^i, are computed from the muon decay rate: 

where |M (p^, p e , pp e , p u ^) \ 2 is the averaged squared amplitude obtained from the Feyn- 
mann diagram at tree level. For polarized muons: 

\M {p^p e ,p Pe ,p^)\ 2 = [fa - m^s)p Pe } {pePu^) , (2) 

where s is the four-spin. For unpolarized muons s — 0. 

g?<3> 3 is the three-body phase-space. In general, the n-body phase space is defined 
by : 

d$ B (P; p x , ...,p n ) = (2tt) 4 5{P -p x - ... - p n ) J] -^-J^ys ■ ( 3 ) 

Differential distributions of decay products are obtained integrating over the phase 
space of the remaining decay particles, 

<rA F (0) (x) + J (0) (x)7Vos# , (4) 



dx d cos 9 



where x denotes the scaled energy, x = 2E e>v / m M and is the average over polarization 
of the initial state muon along the beam direction. 9 is the angle between three- 
momentum of the emitted particle and the muon spin direction and m M is the muon 
mass. The normalized functions and J^°\ in the limit m e = 0, read |16j : 

FW(x) = x 2 (3-2x), Ji 0) (x) =x 2 (l-2x) , (5) 
F<®(x) = x 2 (3 - 2x) , 4°J(x) = x 2 (l - 2x) , (6) 

F£\x) = 6x 2 (l - x) , 4l\x) = 6x 2 (l - x) . (7) 



3 QED corrections 

The QED radiative corrections to the formula (0) where calculated long ago in Ref. [TU] 
through the integration over the neutrino phase-space of the 0(a) corrected differential 
muon decay rate. The QED corrections to the Eq. (jUJ) ( Eq. (J7j) ), are similarly obtained 
from the integration over the v e — e~ [y^ — e~) phase-space. 

In the muon decay, the QED corrected differential rate is given by 

dr = dr + dr y + dr R , (8) 
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where dlY describes the contribution of virtual photon diagrams in Figs. ^3-01 and 
dr^ accounts for the effects of the real photon emission diagrams in Fig. ^ and Fig. 

The virtual photon correction to the decay rate is given by 

dT v = ^ 64G f 2 \M v \ 2 d<5> 3 (p^p e ,p Pe ,pJ , (9) 
where |My| 2 is the squared amplitude. For unpolarized muons it has the expression: 



\M V \ 2 = \M \ 2 - - 
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+ ™ e 9l (P^PPe ) (P^P^ ) + "fy 0R (PePu e ) (p e P^ ) 



(10) 



where g L '^ are ultraviolet (UV) and are listed in the appendix A. The function g[ is 
infrared (IR) divergent while the rest of the u g" functions are finite. 

The IR singularity in g[ is canceled with the soft photon terms of the real emission 
diagrams, which correct the differential rate as 

dT R = -^64G F 2 \M R \ 2 d^{ Pfl ;p e , Pi?e ,p^,k) . (11) 

The explicit expression of the amplitude |M#| 2 can be found in the Appendix A. 

QED corrections for polarized muons are calculated identically to those of unpo- 
larized ones with the replacement p^ — > p^ — sm^ in the amplitudes , where s is the 
muon four- spin [T3] . 



3.1 Soft photon limit and IR cancellation 

Before continuing with the discussion of the exact corrections, let us consider their soft 
photon limit, i.e. k — > 0, as only IR singular terms of virtual and real photon diagrams 
remain in this limit. When soft virtual and soft real photon contributions are added 
up, all 0(a) IR singularities are canceled. 

In this limit, the 0(k) terms in the virtual photon diagrams are neglected and 
Eq. (fTUj) is simplified to : 

\Mr\ 2 =(l-%t) |Mo| 2 , (12) 

where contains all IR divergent terms which are regularized introducing a finite 
photon mass A. 

In the diagrams containing real photon emission, only terms of order 0(k~ 2 ) re- 
main in the soft photon limit. They contain all IR divergent contributions from 
bremsstrahlung. The squared amplitude in Eq.([ll|) reduces to: 

\ M SP ] 2_o 9 a \ _JPj , Pi 2 (PyPe) 1|w|2 / lq \ 
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The divergences in Eq. (112(1 cancel when added with the soft bremsstrahlung part. 
However, Eq. (JT5)l must be previously integrated over the photon-electron phase space 
in order to reduce the real phton emission from a four-body problem to a three-body 
problem. The integral is performed introducing a finite photon mass A, resulting in a 
expression which exactly cancels the IR terms in Eq. (|12|) . 

After the integration over the corresponding phase space, we obtain the soft-photon 
corrected for both and v e distributions, which are proportional to the tree level 
amplitude , in the limits x — > 1 and m e — > 0: 

d 2 N SP rv 

F$ s >(x)[l--£-k(x)]. (14) 



dx d cos 9 v 2n 
The resulting function k(x) is A-independent: 

k(x) = 2L(x) + 27r 2 /3 + ln 2 (l - x). (15) 

where L(x) is the Spence function defined in the Appendix A. 

Realize that k(x) diverges for x — > 1. This singularity is originated due to a failure 
the perturbative treatment: at the end point of the spectrum, the phase space for the 
emission of real photons shrinks to zero and does not compensate the IR infinities of 
the virtual photons. 

The end-point singularity of the corrected electron distribution from muon decay 
has been largely discussed in the literature [TZj. For x — > 1 and k — > 0, the corrected 
electron differential distribution diverges as ln(l — x). Since the IR divergences in 
the muon decay stem from soft-photons in the limit k — > 0, the solution proposed 
to control the end-point divergence is to consider multiple soft-photon emission. The 
effect of considering soft photons at all orders in a is the exponentiation of the singular 
logarithm ln(l — x) which leads to a non-singular distribution [T5] . 

Following as a guideline the solution found for the electron, we apply the same 
proccedure to the neutrino distributions. Consider the neutrino soft-photon correction 
in Eq. (|14|) . At the end-point, for each soft virtual photon and each soft real photon we 
get a ln 2 (l— x) term, which multiplies the tree level amplitude. If there are n soft virtual 
photons and n soft real photons, there are n double logarithms with an additional 
symmetry factor of l/n\. Therefore, the correction to the neutrino distribution at all 
orders in a is obtained summing over n: 

fKX [Fi°\x) + J^V fl (x)cos9]e-^^ 2 ^ . (16) 



dx d cos 9 



The evaluation of infrared divergences at all orders results in the exponentiation of 
the double logarithm, which ensures a non-divergent behavior of the neutrino distribu- 
tions. The exponentiation is only valid for a small region x — > 1. For lower x, we must 
include all the terms of the exact corrections, computed in the next subsection. 



5 



3.2 Results 



Exactly corrected neutrino distributions are obtained considering all terms of the 0(a) 
corrected decay rate, Eq. (JHJ) and integrating over the phase space of the remaining par- 
ticles. Different from the corrected electron distribution, in the neutrino case, the inte- 
grals over the electron-photon phase space in the real emission diagrams are nontrivial. 
We follow the method found in Ref. ^3] to solve analytically these integrals in the cal- 
culation of the QCD corrections to the lepton spectrum from the decay t — > 6 + / + + v\. 
We use the fact that there is a one to one correspondence between the Feynmann dia- 
grams in Fig. for the QED corrections to the //-decay and those for the top quarks. 
This correspondence can be seen by simply replacing 
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a 



-as 



(t, b, /+, pi) 



(17) 



Therefore, by following the techniques detailed Ref. [T4] . we perform the corre- 
sponding phase-space integrals to the differential rate of polarised muons. We find 
that the corrected neutrino angular and energy distributions, including all finite terms 
in the limit m e = 0, are: 



-T^S = F$ (*) + 4? V, (x) cos 9 ^— [F« (x) + JW (x) V, cos 9 
ax a cos 9 ' M 2ir 1 ' ^ 



d 2 N D 



= f£\x) + 4:\x)V,cos9 
ax a cos 9 2tt 



r(o), 



a r 



F£\x) + j£\x)V f ,cos9 



;is) 



where Fp^ v -J^ u are given in Eq. (jHJ), and the one- loop corrections are given by: 



:r 



41' ( 



x 



F£\x) = F^'{x)k{x) + (1 -x) (5 + 8x + 8x 2 ) ln(l - x) 



= F<®(x)k(x) + -(41 - 36x + 42a; 2 - 16x 3 ) ln(l - x) 
» g 

+ -U<82 - 153x + 86x 2 ) , 

= 4° } (x)A;(x) + -(11 - 36x + 14x 2 - 16x 3 - 4/x) ln(l - x) 



+ 



1 

12 1 



-8 + 18x - 103x 2 + 78x d 



(19) 



(20) 



+ -£(10 - 19x) 



(21) 



4l\x) = 4° } (x)A;(2;) + (l-a;)[(-3 + 12a; + 8a; 2 + 4/a;)ln(l-x) 



+ ^ 



2x - 15x 2 ) 



(22) 
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As expected, due to the above correspondence, the results in Eqs. (fT§j) - (|22j) are 
identical to those for the QCD corrections of the lepton distributions from top decay. 

Notice that the function k(x) appears in Eqs. (fT§j) - (|22j) multiplying the tree level 
functions FjP>-jfP\ which agrees with the discussion in the former subsection. 

Fig. |21 and Fig. 0] compare the corrected and the tree level forward and v e distri- 
butions, respectively. In both cases, the relative correction is of 0(0.1%), well below 
the order of the expected precision in the knowledge of the beam parameters. 

The correction found of O(-) ~ 0(0.1%) agrees of that expected from first order 
QED processes. This result differs with the correction of (9(1%) found for the electron 
distribution JU]- The enhacement of the correction the e~ case is due to the "leading 
logs": terms proportional to ln(^ i£ ) which stem from the emission of collinear photons 
in the electron leg JH]- Since neutrinos are not sensitive to QED, no term in -h^^-) 
appears in the neutrino distributions and, neither, terms in -log(^-), which cancel 
when the variables affected by QED corrections, i.e. electron and the photon elec- 
tron momenta, are integrated over. An identical cancellation it is found in the 0(a) 
corrections to the muon lifetime computed which result to be of O(-) ~ (9(0.1%) [TP] . 

In the laboratory frame, neutrino fluxes are boosted along the muon momentum 
direction. The formulae of the corrected distributions in that frame are given in Ap- 
pendix B. 



4 Muon-beam divergence 

We study below the systematic uncertainty in the neutrino distributions produced by 
the muon beam divergence. For the sake of illustration, the quantitative results will be 
given for a 30 GeV unpolarized muon beam decaying in a long straight section pointing 
to a far detector located at 2810 km. 

The natural decay angle of the forward neutrino beam in the laboratory frame is 
deduced from the relation between the rest and laboratory frames. In the rest frame, 
half of the neutrinos are emitted within the cone 9 < tt/2. In the laboratory frame: 

n i cos 9 + 3 . . 

cosd = 1 IT- 23 

l + /3cos6» v ' 



where 3 = y/1 — 7~ 2 is the muon velocity in the laboratory frame. Therefore, half of 
neutrinos are emitted within the cone subtended by the decay angle 9 < I/7. For 
instance, for 30 GeV muons 1/7 = m^/E^ = 3 mrad. 

For the beam and baseline illustrated here, a 10 kt detector and one year of data 
taking [2H], the statistical error on the neutrino flux is of the order of O(0.4%). It 
is then convenient to restrain the uncertainty induced by the muon beam divergence 
below that level. To achieve this, the direction of the beam must be carefully monitored 
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within the decay straight section by placing beam position monitors at its ends. The 
angular divergence of the parent muon beam is then small compared to the natural 
decay angle of the neutrino beam 9 ~ I/7, see Fig. |HJ aiming at present to a divergence 
of 0(0.1/7). It implies that the neutrino beam will be collinear, within the limits set 
by the decay kinematics. 

In our calculations we parameterize this beam focalization by a gaussian distribu- 
tion with standard deviation a ~ O.I/7 (i.e. 0.3 mrad for 30 GeV muon beam) |2*T] . 
which suppresses the flux of neutrinos as they separate from the straight direction. The 
divergence is introduced analytically by considering that the muon direction opens an 
angle a with respect to the z-axis, defined as the direction pointing towards the far 
detector at a distance L, see Fig. El The neutrino distributions in the rest frame, 
Eq. (j3J), are Lorentz boosted along the z-axis. The rest-frame basis (a^cos^) is trans- 
formed to the lab-frame basis (z, cos 8 ), where z = E v / and 9 is the angle between 
the neutrino beam and the z-axis. Using the parameters (3 = y/1 — 7 -2 , the boosted 
distributions read: 

— - v ^ = — ^r— 7 Efz 2 (l — 3(sm (p sin a sin 6 + cos a cos 6 ) ) 
dzdVt nL 2 rrfi, " 

(J, 



x 



I [3m^ — 4zE 2 (l — /3(sin <p sin a sin 6 + cos a cos 6* ))J 
=F Vn \m 2 l — 4z(l — /3(sin<£> sin a sin 6 + cos a cos 6 ))1 j , 



d 2 N - 24-n 
dzdQ 



— —-^t E z (1 — /3(sin (p sin a sin 9 + cos a cos 9 ) ) 
iiL 2 mY, ^ 

x I — 2zE 2 (l — [3(sin.(p sin a sin 9 + cos a cos 9 )) 

T'Pfi m 2 ^ — Az(l — (3(smip sin a sin 9 + cos a cos9 )) j. (24) 
The above expressions are integrated on a, weighted with the gaussian factor 

(25) 



e 2<r^ 



V2na 2 

For unpolarized muons (P M = 0) (for different muon polarizations we obtain similar 
results), it results: 



dzdQ 



An, 



7cL 2 ml " 



AzEl [l-2Pe— cos 6 
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+/3 2 ( 1 e sin 2 Q> sin 2 j + p 2 ( 1 + 6 ~ 2a2 ) cos 2 0' 



d N UeUe - '";/ , • ] _ 2 ) 2 



E*z 2 |mj fl-^e-^cos/) 



-2zE 2 u ( l-2/3e^cos#' 



a* 

?2 1 e I „:„2 /)' • 2 ' , r>2 I 1 + e 



+/? 2 sin 2 fl' sin 2 ^ + ^2 _^ CQS 2 ^ (26) 



Setting 9=0, the expression of forward neutrino fluxes reads: 



dzda tiLV " I » 

'l + e 



-2a 2 ' 



-AzEl 1 - 2/3e^ + 1 



d 2 N UePe 24n M „ 4 2 j ^ 



dzdVt nL 2 m®, M 1 M 



2 



-2(7 2 ' 



-2.Wl-2/^+/3 2 (I±^]H. (27) 



Figs. |H1 show the numerical results for the neutrino and antineutrino spectra in 
a medium baseline (2810 km). We compare the distribution where the muon beam is 
aligned with the detector direction (no beam divergence) with the distribution where 
the muon-beam divergence is included. In the former, neutrino beams are averaged 
over an angle 9 of 0.1 mrad at the far detector [To] . 

Our formulae predict a similar flux correction than previous numerical estimations 
|21j . For instance, a 10% uncertainty in the muon beam divergence would lead to a 
flux uncertainty of 0.3%. We obtain, 

^ Aa 

-Mr ~ °- 03 — 



dE a 



If the muon beam divergence is constrained by lattice design to be less than 0.05/7, 
the loss of flux will be negligible 
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5 Conclusions 



A neutrino factory from muon decay aims at a precision better than (9(1%) in the 
knowledge of the resulting intense neutrino fluxes. 

We have presented here novel results about the effects of QED corrections and muon 
beam divergence on the neutrino differential distributions from muon decay. We have 
given the corresponding corrected formulae (for m e = and m u = 0), including muon 
polarization effects. The induced uncertainties on the neutrino spectra turn out to be 
a safe 0(0.1%). 

Neutrino one-loop corrected distributions diverge at the upper edge of the kinemati- 
cal allowed region. This results from a failure in the cancellation of infrared divergences 
from virtual photons by real photons. Applying the soft photon limit to the exact cal- 
culations, we have isolated the end-point divergent term for the neutrino distributions 
which takes the form of ln 2 (l — x). In order to control this singularity, the double 
logarithmic-contribution is exponentiated, encompassing the contributions from all or- 
ders of perturbation theory. All in all, the exact neutrino distributions get corrections 
of 0(0.1)%, safely below the expected precision in the flux measurements. 

We have also studied carefully the influence of the muon beam divergence on the 
neutrino spectra at the far site. The challenge in designing the neutrino production 
section, where the muons decay, is to constrain the muon beam divergence to a value 
smaller than the natural cone of forward going neutrinos in the laboratory frame, 
(~ I/7). At present, the long straight sections under discussion aim at an angular 
muon beam divergence of the order of O.I/7, typically less than one mrad. 
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Appendix A: QED loop corrections 



A.l Virtual corrections 



There three diagrams containing a photon loop: the exchange of the virtual photon 
between the muon and the electron legs, Fig. [T]d, and lepton propagator corrections, 
Figs. Cfc, ^1. They correct the invariant amplitude of the muon decay as follows: 



-i M 



G 



F r - 



y/2 



{u(p e )T tT u(p v )}{u{p Vll )l' T {l -75MPpJ} • 



(29) 



T a is the corrected \i — e vertex: 

r ff = Tff (i- 75 ) + r»+r^, (30) 

where results from the diagram in Fig. [I]p and T c j d from those in Figs, ^p, ^i. 

After integration over the photon momentum, the correction from diagram [I]p, has 
the expression 



a 



-7T [ km + <7uv) 7. (1 - 7s) + 9r la (1 + 75) 
Zix 



+ ^lPi ct (1-75)+^rP2 CT (1 + 7s) 



whith 



9m 



coth< 



L 



L 



2 sinh ( 



e<P - e" 



1 sinh 1 



2 (cosh u; — cosh 



2 sinh 



e w — e - 



cusinhcu 1 / , / A 2 
+ - w - In — 
2 I Im 2 . 



In 



'2 sinh (^) 
,2 sinh (*±4 



3 
2' 



sinh 1 



2m M sinh < 



2m P sinh i 



sinh 1 



1 sinh 00 



+ 



cosh to — cosh 
ui sinh <f) — <p sinh a; 
cosh u; — cosh 



(31) 



(32) 



where the IR term is regularized by a finite photon mass A and the variables 

(P/*Pe) u _ TO,, 



cosh( 



TO M m e 



m e 



(33) 



are introduced following [TU] . 
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L(x) is the Spence function 

L(x) 



dt 



In |1 -t\ 

t 



(34) 



The contribution of self-energy diagrams to the muon-electron vertex, after inte- 
gration over the photon momentum, is given by 

a 1 



( h xjv + h m) 



2tt 2 



where, now, 



h 



c.d 
UV 



-i L-ln £ 



h 



c.d 
IK 



2 

uj — In 



rat 



A 2 



rat 



2 . 



(35) 

(36) 
(37) 



Adding Eq. (J3~2*j) and Eq. (jHSj) the UV divergences are exactly cancelled. The IR 
terms in Eq. (|37j) . when combined with Eq. (|3*T|) gives rise to the term 



<7l 



+ 



• coth i 
i sinh a 



L 



2 sinh < 



L 



2 sinh < 



ye" 7 — e~ 
a; sinh 



2 (cosh a; — cosh( 



+ 2-1 



• coth < 



uj — In 



(uj — 0) In 
/A 2 \ \ 



'2 sinh 
2 sinh 



(3* 



A. 2 Bremstrahlung corrections 

The contribution from real photon emission, Fig. ^ and Fig. [If, is given by 

1 



2m, 



64 Gf |M fl | d$4(pfi;p e ,pp e ,p v , k), 



where the amplitude |Mr| 2 has the following expression: 



Mr\ = r— 



a 
2^ 



A 



+ 



B 



C 



.(Pf,k) 2 (p e k) 2 {p^k)(p e k)_ 
The numerators for unpolarized muons read: 

A = Pi [{P^Pve ) (PePuu ) - {kpp e ) (pePup ) - (jPyk) {kp Pe ) {PeP^)} 



B 

C 



(39) 



(40) 



Pi [(PvPPe)(PePv„) + (P»PPe)(kpvJ - (P»Pi? e )(Pek)(kp Vfl )] 
{P^Pe) [ZiPuPvMPePvJ + (PvPP e ){kPu„) ~ {kpp e ){PePu^)] (41) 
+ (Pek) [{PfiPv e ) ) + (p^ ) ) - (p M Pp e ) (&Pi/„)] 

- (P^KPf.PuJiPePu^) + {PePu e ){PeP^) ~ (fcp P J (PeP^)] ■ 

Terms of order fc 2 are not included in Eq. (|41j) , since they vanish in the limit of massless 
photons. 
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Appendix B: QED corrected distributions in the 
laboratory frame 

In order to obtain the neutrino distributions in the laboratory frame, a Lorentz 
boost is performed in the direction of the muon velocity towards the detector at a dis- 
tance L. The rest-frame basis (x, cos 9) is transformed to the lab- frame basis (z, cos 9 ), 
where z = E u /E^ is the scaled energy at the laboratory frame and 9 is the angle be- 
tween the neutrino beam and the direction of the muon beam ^3]. The muon beam 
divergence is set to zero. Using the parameters 7 = E^/m^ and (3 = y/1 — 7~ 2 , the 
boosted distributions read: 



where 

4%, 9' 



d 2 N v 



dzd cos 9' 



F^(z,9')+V,4°J(z,9')cos9' 
F^(z,9')+V,4 1 J(z 1 9')cos9' 



a 
2^ 



d 2 N D 



F£\z,9')+V,4: > (z,9') C o S 9 



dzd cos 9' 



(0), 



a 
2^ 



F£\z,9')+V, (z, 9') cos 9 



E 



-^z 2 (l-f3cos9')(3m 2 -4Elz(l-f3cos9')), 



rrr, 



E 



8 -f z 2 ( 1 - /3 cos 9' ) (m 2 - AE' 2 z( 1 - {3 cos 0' ) ) , 



??? 



i? 4 

48^- z 2 (1 - pcos9')(ml - 2E 2 ^(1 -0008 0% 
£ 4 

48-f z 2 (1 - cos 6') (ml - 2£ 2 z(l - Pcos0 )), 



F^(z,9')k(z,9') + 



3(1 -p cos 0') 
+ 42 (2 7 2 z(l -/Jcosfl')) 2 - 16 (2 7 2 z(l - pcos9') 

x ln(l - 2 7 2 z(l -/3cos0')) 

+ - (2 7 2 z(l -/3cos0')) [82 - 153 (2 7 2 z(l-/3cos^' 
+ 86 (2 7 2 z(l -/^costf' 



(42) 

(43) 

(44) 
(45) 
(46) 
(47) 



41 - 36 (2Yz (1-/3 cos 9 
3" 



(4* 
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4?W) = 4?M>M') + 



11 - 36 (27 2 2(1 -Pcos 6') 



f£\z,9') = F£\ z ,e')k(z,e') + 



3(1 — (3 cos 9') 
+ 14 (2-f 2 z {1-13 cos e')f - 16 fa 2 z (1 - (3 cos 6 ' )f 

+ A(2-i 2 z(l- (3cos6')y l x ln(l - 2 7 2 ^(1 -/3cosfl 
8 + 18 (2 7 2 ^(1 -/3cos0')) 
103 (2 7 2 ^ (1 - /^costf')) 2 + 78 (2 7 2 z (1 - /3cos#')) 3 

2 (l -2 7 2 ^(1 -pcosO' 



+ 



1 r 
2 



X 



X 



(l-/3cos0') 
5 + 8 (2 7 2 ,z(l-/3cos0')) 

+ 8 (2 7 2 ^(1 -{3cos6')f ln(l - 2 7 2 ^ (1 -/3cos0')) 

+ ^ (2 7 2 z(l -/3cos0')) [lO- 19 (2 7 2 z(l-/3cos0' 



, 2(l- 2 7 2 z(l-/3cos0' 



(1 -/3cos0') 

-3 + 12 (2 7 2 ^ (1 - /3 cos 0'j) + 8 (2 7 2 ^ (1 - /?cos0'))' 
+ 4 (2 7 2 z(l-/?cos0')) _1 ln(l - 2 7 2 z(l -/3cos0') 



+ - 



8 - 2 (2 7 2 ^ (1-/3 cos 0')) - 15 (2 7 2 ^ (1-/3 cos 0' 



(49) 



(50) 



(51) 



k(z,9') = ln 2 (l- 2 7 2 ^(l-/3cos0')) + 2L(2 7 2 z (1-/3 cos 0' 



2tt 2 



(52) 
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Figure 1: QED Radiative corrections to muon decay: (a) Tree level, (b) Vertex cor- 
rection, (c) Muon propagator correction, (d) Elctron propagator correction, (e) Muon leg 
bremsstrahlung, (d) Elctron leg bremsstrahlung 
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;ure 2: Zeroth order and 0(a) corrected forward distributions for = 30 GeV and 
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Figure 3: Detail of Fig. El 
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Figure 5: Detail of Fig. H 
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Figure 7: and differential distributions. The solid lines represent the spectra obtained by 
averaging over an angular divergence of 0.1 mrad and the dashed lines the spectra including muon 
beam divergence. The distributions are plotted in the forward direction cos# = pointing towards 
a detector located 2810 km from a the neutrino source of unpolarized positive or negative muons 
circulating in the storage ring with energies of 30 GeV. 
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Fi gure 8: v e and v e differential distributions. The solid lines represent the spectra obtained by 
averaging over an angular divergence of 0.1 mrad and the dashed lines the spectra including muon 
beam divergence. The distributions are plotted with the same parameters as of figi[7|l. 
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